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ABSTRACT.  During  the  past  five  to  ten  years,  a  broad  partnership  of  institutions 
under  NOPP  sponsorship  has  collaborated  in  developing  and  demonstrating  the 
performance  and  application  of  eddy- resolving,  real-time  global-  and  basin-scale 
ocean  prediction  systems  using  the  HYbrid  Coordinate  Ocean  Model  (HYCOM). 

The  partnership  represents  a  broad  spectrum  of  the  oceanographic  community, 
bringing  together  academia,  federal  agencies,  and  industry/commercial  entities,  and 
spanning  modeling,  data  assimilation,  data  management  and  serving,  observational 
capabilities,  and  application  of  HYCOM  prediction  system  outputs.  In  addition  to 
providing  real-time,  eddy-resolving  global-  and  basin-scale  ocean  prediction  systems 
for  the  US  Navy  and  NOAA,  this  project  also  offered  an  outstanding  opportunity  for 
NOAA-Navy  collaboration  and  cooperation,  ranging  from  research  to  the  operational 
level.  This  paper  provides  an  overview  of  the  global  HYCOM  ocean  prediction  system 
and  highlights  some  of  its  achievements.  An  important  outcome  of  this  effort  is 
the  capability  of  the  global  system  to  provide  boundary  conditions  to  even  higher 
resolution  regional  and  coastal  models. 


INTRODUCTION 

A  broad  partnership  of  institutions  has 
collaborated  over  the  past  five  to  ten 
years  to  develop  and  demonstrate  the 
performance  and  application  of  eddy¬ 
resolving,  real-time  global-  and  basin 
scale  ocean  prediction  systems  using 
the  HYbrid  Coordinate  Ocean  Model 
(HYCOM).  These  systems  are  in  the 
process  of  being  transitioned  to  opera¬ 
tional  use  by  the  US  Navy  at  the  Naval 
Oceanographic  Office  (NAVOCEANO), 
Stennis  Space  Center,  MS,  and  by  the 
National  Oceanic  and  Atmospheric 
Administration  (NOAA)  at  the  National 
Centers  for  Environmental  Prediction 
(NCEP),  Washington,  DC.  The  systems 
run  efficiently  on  a  variety  of  massively 
parallel  computers  and  include  sophis 
ticated,  but  relatively  inexpensive,  data 
assimilation  techniques  for  satellite 
altimeter  sea  surface  height  and  sea 
surface  temperature  as  well  as  in  situ 
temperature,  salinity,  and  float  displace¬ 
ment.  The  partnership  represents  a  broad 
spectrum  of  the  oceanographic  commu¬ 
nity,  bringing  together  academia,  federal 


agencies,  and  industry/commercial 
entities  in  activities  that  span  modeling, 
data  assimilation,  data  management 
and  serving,  observational  capabilities, 
and  application  of  HYCOM  predic¬ 
tion  system  outputs.  All  participating 
institutions  were  committed  and  the 
collaborative  partnership  provided  an 
opportunity  to  leverage  and  accelerate 
the  efforts  of  existing  and  planned 
projects,  consequently  producing  a  high- 
quality  product  that  should  collectively 
serve  a  wider  range  of  users  than  would 
the  individual  projects. 

The  collaboration  was  initiated  in  the 
late  1990s  by  ocean  modelers  at  the  Naval 
Research  Laboratory,  Stennis,  MS,  who 
approached  colleagues  at  the  University 
of  Miami’s  Rosenstiel  School  of  Marine 
and  Atmospheric  Science  regarding  an 
extension  of  the  range  of  applicability  of 
the  US  Navy  operational  ocean  predic¬ 
tion  system  to  coastal  regions  (e.g.,  the 
US  Navy  systems  at  the  time  were  seri¬ 
ously  limited  in  shallow  water  and  in 
handling  the  transition  from  deep  to 
shallow  water).  HYCOM  (Bleck,  2002) 


was  therefore  designed  to  extend  the 
range  of  existing  operational  Ocean 
General  Circulation  Models  (OGCMs). 
The  freedom  to  adjust  the  vertical  spacing 
of  the  generalized  (or  hybrid)  coordinate 
layers  in  HYCOM  simplifies  the  numer 
leal  implementation  of  several  processes 
and  allows  for  a  smooth  transition 
from  the  deep  ocean  to  coastal  regimes. 
HYCOM  retains  many  of  the  charac¬ 
teristics  of  its  predecessor,  the  Miami 
Isopycnic  Coordinate  Ocean  Model 
(Bleck  et  al,  1992;  Bleck  and  Chassignet, 
1994),  while  allowing  coordinates  to 
locally  deviate  from  isopycnals  wherever 
the  latter  may  fold,  outcrop,  or  generally 
provide  inadequate  vertical  resolution. 
The  collaboration  led  to  the  development 
of  a  consortium  for  hybrid-coordinate 
data  assimilative  ocean  modeling 
supported  by  NOPP  to  make  HYCOM 
a  state-of-the-art  community  ocean 
model  with  data  assimilation  capability 
that  could:  (1)  be  used  in  a  wide  range  of 
ocean-related  research,  (2)  become  the 
next  generation  eddy-resolving  global 
ocean  prediction  system,  and  (3)  be 
coupled  to  a  variety  of  other  models, 
including  littoral,  atmospheric,  ice,  and 
biochemical  models.  One  outcome  of  this 
collaboration  was  the  establishment  of  a 
near-real-time  North  Atlantic  prediction 
system  based  on  HYCOM. 

In  2003,  NOAA  NCEP  joined  forces 
with  the  HYCOM  consortium  in 
responding  to  a  NOPP  Broad  Agency 
Announcement  aimed  at  “imple¬ 
menting  the  initial,  preoperational 
US  contribution(s)  to  the  Global 
Ocean  Data  Assimilation  Experiment 
(GODAE).”  These  efforts  were  intended 
to  be  pilot  projects  under  Ocean. US, 
the  National  Office  for  Integrated  and 
Sustained  Ocean  Observations,  and  to 
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eventually  lead  to  sustained  operational 
efforts  supported  by  US  agencies  such 
as  NOAA  and  the  Navy.  The  HYCOM 
consortium  therefore  became  one  of  the 
US  components  of  GODAE,  a  coordi¬ 
nated  international  system  of  observa¬ 
tions,  communications,  modeling, 
and  assimilation  that  delivers  regular, 
comprehensive  information  on  the 
state  of  the  ocean  (see  Chassignet  and 
Verron,  2006,  for  a  review).  Navy  and 
NOAA  applications,  such  as  maritime 


safety,  fisheries,  the  offshore  industry, 
and  management  of  shelf/coastal  areas, 
are  among  the  expected  beneficiaries  of 
the  HYCOM  ocean  prediction  systems 
(http://www.hycom.org).  More  specifi¬ 
cally,  the  precise  knowledge  and  predic¬ 
tion  of  ocean  mesoscale  features  helps 
the  Navy,  NOAA,  the  Coast  Guard, 
industry,  and  fisheries  with  endeavors 
such  as  ship  and  submarine  routing, 
search  and  rescue,  oil  spill  drift  predic¬ 
tion,  open-ocean  ecosystem  monitoring, 


fisheries  management,  and  short-range 
coupled  atmosphere- ocean,  coastal,  and 
nearshore  environmental  forecasting. 

In  addition  to  operational  eddy¬ 
resolving  global-  and  basin-scale  ocean 
prediction  systems  for  the  US  Navy 
and  NOAA,  respectively,  this  project 
offered  an  outstanding  opportunity  for 
NOAA-Navy  collaboration  and  coop¬ 
eration  ranging  from  research  to  the 
operational  level. 

BACKGROUND 

Over  the  past  several  decades,  numerical 
modeling  studies  have  demonstrated 
progress  in  both  model  architecture 
and  the  availability  of  computational 
resources  to  the  scientific  community. 
Perhaps  the  most  noticeable  aspect  of 
these  advances  has  been  the  evolution 
from  simulations  on  coarse-resolution 
horizontal/vertical  grids  outlining  basins 
of  simplified  geometry  and  bathymetry 
and  forced  by  idealized  stresses,  to  fine- 
resolution  simulations  incorporating 
realistic  coastal  definition  and  bottom 
topography  and  forced  by  observational 
data  on  relatively  short  time  scales 
(Hurlburt  and  Hogan,  2000;  Smith  et  al., 
2000;  Chassignet  and  Garraffo,  2001; 
Maltrud  and  McClean,  2005;  Hurlburt 
et  al.,  2008).  The  choice  of  the  vertical 
coordinate  system,  however,  remains 
one  of  the  most  important  aspects  of  an 
ocean  models  design.  In  practice,  the 
representation  and  parameterization  of 
processes  not  resolved  by  the  model  grid 
are  often  directly  linked  to  the  vertical 
coordinate  choice  (Griffies  et  al.,  2000). 
Oceanic  general  circulation  models 
traditionally  represent  the  vertical  in 
a  series  of  discrete  intervals  in  either  a 
depth,  density,  or  terrain- following  unit. 
Because  none  of  the  three  main  vertical 
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coordinates  (depth,  density,  and  terrain¬ 
following)  provide  universal  optimality,  it 
is  natural  to  envision  a  hybrid  approach 
that  combines  the  best  features  of  each 
vertical  coordinate.  Isopycnic  (potential 
density- tracking)  layers  work  best  for 
modeling  the  deep  stratified  ocean;  levels 
at  constant  fixed  depth  or  pressure  are 
best  for  providing  high  vertical  resolution 
near  the  surface  within  the  mixed  layer; 
and  terrain-following  levels  are  often  the 
best  choice  for  modeling  shallow  coastal 
regions.  In  HYCOM,  the  optimal  vertical 
coordinate  distribution  of  the  three 
vertical  coordinate  types  is  chosen  at 
every  time  step  and  in  every  grid  column 
individually.  The  default  configuration  of 
HYCOM  is  isopycnic  in  the  open  strati¬ 
fied  ocean,  but  it  makes  a  dynamically 
and  geometrically  smooth  transition  to 
terrain -following  coordinates  in  shallow 
coastal  regions  and  to  fixed  pressure -level 
coordinates  in  the  surface  mixed  layer 
and/or  unstratified  open  seas.  In  doing 
so,  the  model  takes  advantage  of  the 
different  coordinate  types  in  optimally 
simulating  coastal  and  open-ocean  circu¬ 
lation  features  (Chassignet  et  al.,  2003, 
2006,  2007).  A  user-chosen  option  allows 
specification  of  the  vertical  coordinate 
separation  that  controls  the  transition 
among  the  three  coordinate  systems. 

The  assignment  of  additional  coordinate 
surfaces  to  the  oceanic  mixed  layer  also 
allows  the  straightforward  implementa¬ 
tion  of  multiple  vertical  mixing  turbu¬ 
lence  closure  schemes  (Halliwell,  2004). 
The  choice  of  the  vertical  mixing  param  ¬ 
eterization  is  also  of  importance  in  areas 
of  strong  entrainment,  such  as  overflows. 

Data  assimilation  is  essential  for 
ocean  prediction  because:  (a)  many 
ocean  phenomena  are  due  to  nonlinear 
processes  (i.e.,  flow  instabilities)  and 


thus  are  not  a  deterministic  response  to 
atmospheric  forcing,  (b)  errors  exist  in 
the  atmospheric  forcing,  and  (c)  ocean 
models  are  imperfect,  including  limita¬ 
tions  in  numerical  algorithms  and  in 
resolution.  Most  of  the  information  about 
the  ocean  surfaces  space-time  variability 
is  obtained  remotely  from  instruments 
aboard  satellites  (sea  surface  height  [SSH] 
and  sea  surface  temperature  [SST]),  but 
these  observations  are  insufficient  for 
specifying  the  subsurface  variability. 
Vertical  profiles  from  expendable 
bathythermographs  (XBT),  conductivity- 
temperature  depth  (CTD)  profilers, 
and  profiling  floats  (e.g.,  Argo,  which 
measures  temperature  and  salinity  in 
the  upper  2000  m  of  the  ocean)  provide 
another  substantial  source  of  data.  Even 
together,  these  data  sets  are  insufficient 
to  determine  the  state  of  the  ocean 
completely,  so  it  is  necessary  to  exploit 
prior  statistical  knowledge  based  on 
past  observations  as  well  as  our  present 
understanding  of  ocean  dynamics.  By 
combining  all  of  these  observations 
through  data  assimilation  into  an  ocean 
model,  it  is  possible,  in  principle,  to 
produce  a  dynamically  consistent  depic¬ 
tion  of  the  ocean.  However,  in  order 
to  have  any  predictive  capabilities,  it 
is  extremely  important  that  the  freely 
evolving  ocean  model  (i.e.,  non-data 
assimilative  model)  is  skilled  in  repre¬ 
senting  ocean  features  of  interest. 

To  properly  assimilate  the  SSH  anoma¬ 
lies  determined  from  satellite  altimeter 
data,  the  oceanic  mean  SSH  over  the 
altimeter  observation  period  must  be 
provided.  In  this  mean,  it  is  essential  that 
the  mean  current  systems  and  associated 
SSH  fronts  be  accurately  represented  in 
terms  of  position,  amplitude,  and  sharp¬ 
ness.  Unfortunately,  Earths  geoid  is  not 


presently  known  with  sufficient  accuracy 
for  this  purpose,  and  coarse  hydrographic 
climatologies  (~  0.5°-l°  horizontal 
resolution)  cannot  provide  the  spatial 
resolution  necessary  when  assimilating 
SSH  in  an  eddy- resolving  model  (hori  ¬ 
zontal  grid  spacing  of  1/10°  or  finer).  At 
these  scales  of  interest,  it  is  essential  to 
have  the  observed  means  of  boundary 
currents  and  associated  fronts  sharply 
defined  (Hurlburt  et  al.,  2008).  Figure  1 
shows  the  climatological  mean  derived 
on  a  0.5°  grid  using  surface  drifters  by 
Maximenko  and  Niiler  (2005)  as  well 
as  the  mean  currently  used  in  the  Navy 
global  HYCOM  prediction  system 
(see  following  section  for  details).  The 
HYCOM  mean  was  constructed  as 
follows:  a  five-year  mean  SSH  field  from 
a  non-data-assimilative  1/12°  global 
HYCOM  run  was  compared  to  available 
climatologies,  and  a  rubber-sheeting 
technique  (Carnes  et  al.,  1996)  was  used 
to  modify  the  model  mean  in  two  regions 
(the  Gulf  Stream  and  the  Kuroshio) 
where  the  western  boundary  current 
extensions  were  not  well  represented 
and  where  an  accurate  frontal  location 
is  crucial  for  ocean  prediction.  Rubber 
sheeting  involves  a  suite  of  computer 
programs  that  operate  on  SSH  fields, 
overlaying  contours  from  a  reference  field 
and  moving  masses. 

THE  HYCOM  OCEAN 
PREDICTION  SYSTEMS 

Two  systems  are  currently  run  in  real 
time  by  the  US  Navy  at  NAVOCEANO, 
Stennis  Space  Center,  MS,  and  by 
NOAA  at  NCEP,  Washington,  DC 
(http://www.hycom.org). 

The  first  system  is  the  NOAA  Real 
Time  Ocean  Forecast  System  for  the 
Atlantic  (RTOFS- Atlantic),  which  has 
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Figure  1.  The  top  panel  shows  mean  sea  surface  height  (SSH;  in  cm)  derived  from  surface  drifters 
(Maximenko  and  Niiler,  200S),  and  the  bottom  panel  shows  the  same  from  a  non  data  assimila 
five  HYbrid  Coordinate  Ocean  Model  (HYCOM)  run  corrected  in  the  Gulf  Stream  and  Kuroshio 
regions  using  a  rubber  sheeting  technique.  The  RMS  difference  between  the  two  fields  is  9  2  cm 


been  running  in  real  time  since  2005. 

The  Atlantic  domain  spans  25°S  to  76°N 
with  a  horizontal  resolution  varying 
from  4  km  near  the  US  coastline  to 
20  km  near  the  African  coast.  The  system 
is  run  daily  with  one-day  nowcasts 
and  five-day  forecasts.  Prior  to  June 
2007,  only  the  sea  surface  temperature 
was  assimilated.  In  June  2007,  NOAA 
implemented  the  3D-variational  data 
assimilation  of:  (1)  sea  surface  tempera¬ 
ture  and  sea  surface  height  (Jason- 1, 
Geosat  Follow-On  [GFO],  and  soon 
Envisat),  (2)  temperature  and  salinity 
profile  assimilation  (e.g.,  Argo,  CTDs, 
moorings),  and  (3)  GOES  data.  Plans 


are  to  expand  this  system  globally  using 
the  US  Navy  configuration  described 
in  the  following  paragraph.  The  NCEP 
RTOFS-Atlantic  model  data  is  distrib¬ 
uted  in  real  time  through  NCEPs 
operational  ftp  server  (ftp://ftpprd.ncep. 
noaa.gov)  and  the  NOAA  Operational 
Model  Archive  and  Distribution  System 
(NOMADS;  http://nomads6.ncdc.noaa. 
gov/ncep_data/index.html)  server.  The 
latter  server  is  also  using  Open  Project 
for  a  Network  Data  Access  Protocol 
(OPeNDAP)  middleware  as  a  data-access 
method.  NCEPs  RTOFS  Atlantic  model 
data  is  also  archived  at  the  National 
Oceanographic  Data  Center  (NODC, 


http://data.nodc.noaa.gov/ncep/rtofs). 

The  second  system  is  the  global  US 
Navy  nowcast/forecast  system  using 
the  1/12°  global  HYCOM  (6.5  km  grid 
spacing  on  average,  3.5-km  grid  spacing 
at  the  North  Pole,  and  32  hybrid  layers 
in  the  vertical),  which  has  been  running 
in  near  real  time  since  December  2006 
and  in  real  time  since  February  2007. 

The  current  ice  model  is  thermody¬ 
namic,  but  it  will  soon  include  more 
physics  as  it  is  upgraded  to  the  Polar  Ice 
Prediction  System  (PIPS,  based  on  the 
Los  Alamos  ice  model  known  as  CICE). 
The  model  is  currently  running  daily  on 
379  processors  on  an  IBM  Power  5+  at 
NAVOCEANO  using  a  part  of  the  opera¬ 
tional  allocation  on  the  machine.  The 
daily  run  consists  of  a  five-day  hindcast 
and  a  five-day  forecast  and  takes  about  1 5 
wall  clock  hours.  The  system  assimilates 

(1)  SSH  (Envisat,  GFO,  and  Jason- 1), 

(2)  SST  (all  available  satellite  and  in 
situ  sources),  (3)  all  available  in  situ 
temperature  and  salinity  profiles  (e.g., 
Argo,  CTDs,  moorings),  and  (4)  Special 
Sensor  Microwave/ Imager  (SSMI)  sea  ice 
concentration.  The  three-dimensional 
multivariate  optimum  interpolation 
Navy  Coupled  Ocean  Data  Assimilation 
(NCODA)  system  (Cummings,  2005)  is 
the  assimilation  technique.  The  NCODA 
horizontal  correlations  are  multivariate 
in  geopotential  and  velocity,  thereby 
permitting  adjustments  (increments) 

to  the  mass  field  to  be  correlated  with 
adjustments  to  the  flow  field.  The 
velocity  adjustments  are  in  geostrophic 
balance  with  the  geopotential  incre¬ 
ments,  and  the  geopotential  increments 
are  in  hydrostatic  agreement  with  the 
temperature  and  salinity  increments. 
Either  the  Cooper  and  Haines  (1996) 
technique  or  synthetic  temperature  and 
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Figure  2.  Verification  of  30  day  ocean  forecasts:  median  SSH  anomaly  correlation  vs.  forecast  length  in  comparison  with  the  verifying  analysis  for  the 
global  US  Navy  HYCOM  over  the  world  ocean  and  five  subregions.  The  red  curves  verify  forecasts  using  operational  atmospheric  forcing,  which  reverts 
toward  climatology  after  five  days.  The  green  curves  verify  "forecasts"  with  analysis  quality  forcing  for  the  duration,  and  the  blue  curves  verify  forecasts 
of  persistence  (i.e.;  no  change  from  the  initial  state).  The  plots  show  median  statistics  over  twenty  30  day  HYCOM  forecasts  initialized  during  January 
2004 -December  2005,  a  period  when  data  from  three  nadir-beam  altimeters,  Envisat,  GEOSAT  Follow-on,  and  Jason- 1,  were  assimilated.  The  reader  is 
referred  to  Hurlburt  et  al.  (2008)  and  an  article  scheduled  for  the  September  2009  issue  of  Oceanography  for  a  more  detailed  discussion  of  these  results. 


salinity  profiles  (Fox  et  al.,  2002)  can  be 
used  for  downward  projection  of  SSH 
and  SST.  Figure  2  shows  an  example  of 
forecast  performance. 

Validation  of  the  results  is  underway 
using  independent  data  with  a  focus  on 
the  large-scale  circulation  features,  SSH 
variability,  eddy  kinetic  energy,  mixed- 
layer  depth,  vertical  profiles  of  tempera¬ 
ture  and  salinity,  SST,  and  coastal  sea 
levels.  Figures  3  and  4  show  examples  for 
the  Gulf  Stream  region  while  Figure  5 
documents  the  performance  of  HYCOM 
in  representing  the  mixed-layer  depth. 
HYCOM  is  also  an  active  participant  in 
the  international  GODAE  comparison  of 
global  ocean  forecasting  systems. 


DISTRIBUTION  OF 
GLOBAL  HYCOM  HINDCASTS 
AND  FORECASTS 

The  model  outputs  from  the  global 
US  Navy  hindcast  experiment  from 
November  2003  to  present  are  avail¬ 
able  through  the  HYCOM  consortium 
Web  page,  http://www.hycom.org.  The 
HYCOM  data  distribution  team  devel¬ 
oped  and  implemented  a  comprehensive 
data  management  and  distribution 
strategy  that  allowed  easy  and  efficient 
access  to  the  global  HYCOM-based 
ocean  prediction  system  output  to 
(a)  coastal  and  regional  modeling 
groups,  (b)  the  wider  oceanographic  and 
scientific  community,  including  climate 


and  ecosystem  researchers,  and  (c)  the 
general  public.  The  outreach  system 
consists  of  a  Web  server  that  acts  as  a 
gateway  to  backend  data  management, 
distribution,  and  visualization  applica¬ 
tions  (http://www.hycom.org/dataserver). 
These  applications  enable  end  users  to 
obtain  a  broad  range  of  services  such  as 
browsing  of,  for  example,  data  sets,  GIF 
images,  NetCDF  files,  and  FTP  requests 
of  data.  The  100  terabytes  HYCOM  Data 
Sharing  System  is  built  upon  two  existing 
software  components:  OPeNDAP  (see 
article  by  Cornillon  et  al.,  this  issue) 
and  the  Live  Access  Server  (LAS;  http:// 
ferret.pmel.noaa.gov/LAS/).  These  tools 
and  their  data  distribution  methods  are 
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Figure  3.  Surface  (top  panels)  and  700  m  (lower  panels)  eddy  kinetic  energy  from  observations  (left  panels)  and  HYCOM  (right 
panels)  over  the  period  2004-2006.  The  observed  surface  eddy  kinetic  energy  (upper  left  panel)  is  from  Fratantoni  (2001)  and  the 
700  m  eddy  kinetic  energy  (lower  left  panel)  is  from  Schmitz  (1996).  The  units  are  in  cm  s  The  Gulf  Stream  north  wall  position 
±  1  standard  deviation  is  overlaid  on  the  top  panels. 
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Figure  4.  Modeled  analysis  of  the  sea 
surface  height  field  on  September  8,  2008. 
The  white  line  represents  the  independent 
frontal  analysis  of  sea  surface  temperature 
observations  performed  by  the  Naval 
Oceanographic  Office. 
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Figure  5  Median  bias  error  (in  m)  of  mixed  layer  depth  (MLD)  calculated  from  simulated  and  approximately  66,000  unassimilated 
observed  profiles  over  the  period  June  2007-May  2008.  Blue  (red)  indicates  a  simulated  MLD  shallower  (deeper)  than  observed; 
S3%  of  the  simulated  MLDs  are  within  10  m  of  the  observations,  and  these  are  represented  as  gray.  The  basinwide  median  bias 
error  is  6.6  m  and  the  RMS  error  is  40  m. 


described  below.  In  the  current  setup, 
the  OPeNDAP  component  provides  the 
middleware  necessary  to  access  distrib¬ 
uted  data,  while  the  LAS  functions  as  a 
user  interface  and  a  product  server.  The 
abstraction  offered  by  the  OPeNDAP 
server  also  makes  it  possible  to  define  a 
virtual  data  set  that  LAS  will  act  upon, 
rather  than  physical  files.  An  OPeNDAP 
“aggregation  server”  uses  this  approach 
to  append  model  time  steps  from  many 
separate  files  into  virtual  data  sets.  The 
HYCOM  Data  Service  has  been  in  opera¬ 
tion  for  the  last  four  years  and  has  seen 
a  steady  increase  in  the  user  base.  In  the 
last  year,  the  service  received  approxi¬ 
mately  20,000  hits  per  month.  In  addition 
to  the  numerous  requests  from  educa¬ 
tional  institutions  and  researchers,  this 
service  has  been  providing  near-real-time 


data  products  to  several  private  compa 
nies  in  France,  the  Netherlands,  Portugal, 
and  the  United  States. 

BOUNDARY  CONDITIONS  FOR 

REGIONAL  AND  COASTAL 
MODELS  NESTED  IN  HYCOM 

An  important  attribute  of  the  data 
assimilative  HYCOM  system  is  its  capa 
bility  to  provide  boundary  conditions 
to  even  higher-resolution  regional  and 
coastal  models.  The  current  horizontal 
and  vertical  resolution  of  the  global 
forecasting  system  marginally  resolves 
the  coastal  ocean  (7  km  at  mid  latitudes, 
with  up  to  15  terrain- following  coordi 
nates  over  the  shelf)*  but  it  is  an  excellent 
starting  point  for  even  higher-resolution 
coastal  ocean  prediction  efforts.  Several 
partners  within  the  HYCOM  consortium 


evaluated  the  boundary  conditions  and 
demonstrated  the  value  added  by  the 
global  and  basin  HYCOM  data  assimila¬ 
tive  system  output  for  coastal  ocean 
prediction  models.  The  inner  nested 
models  may  or  may  not  be  HYCOM 
(i.e.,  the  nesting  procedure  can  handle 
any  vertical  grid  choice).  Outer  model 
fields  are  interpolated  to  the  horizontal 
and  vertical  grid  of  the  nested  model 
throughout  the  entire  time  interval  of 
the  nested  model  simulation  at  a  time 
interval  specified  by  the  user,  typically 
once  per  day.  The  nested  model  is  initial 
lzed  from  the  first  archive  file  and  the 
entire  set  of  archives  provides  boundary 
conditions  during  the  nested  run, 
insuring  consistency  between  initial  and 
boundary  conditions.  This  procedure  has 
proven  to  be  very  robust. 
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Figure  6.  Sea  surface  temperature  (fiC)  and  surface  velocity  fields  from  the  Regional 
Ocean  Modeling  System  (ROMS)  West  Florida  Shelf  domain  (inside  the  dashed 
lines)  and  the  F4YCOM  ocean  prediction  system  (outside  the  dashed  lines). 


Wf-S  nested  ri  Climatology  (111  exp  )  WFS  nested  in  Hycom  T  and  S  {2nd  exp  ) 


Figure  7.  Mean  sea  surface  height  (in  in)  on  April  2004  from  the  model  (first  three  panels)  and 
observations  (lower  right  panel). 


We  now  describe  examples  that  illus¬ 
trate  applications  that  benefited  from  the 
successful  partnership  between  coastal 
ocean  modelers  and  ocean  prediction 
system  developers.  The  main  goal  of  the 
first  coastal  application  was  to  assess  the 
impact  of  open  boundary  conditions 
from  the  HYCOM  ocean  prediction 
system  on  the  dynamics  and  accuracy 
of  a  regional  West  Florida  Shelf  (WFS) 
model  (an  effort  led  by  Robert  Weisherg 
at  the  University  of  South  Florida).  The 
Gulf  of  Mexico  Loop  Current  is  the  main 
large-scale  ocean  feature  in  the  WFS 
domain  as  illustrated  by  Figure  6,  which 
shows  an  example  of  the  SST  and  surface 
velocity  fields  from  the  West  Florida 
Shelf  domain  embedded  in  the  US  Navy 
HYCOM  ocean  prediction  system  (Barth 
et  al.,  2008).  The  ocean  model  consists 
of  the  Regional  Ocean  Modeling  System 
(ROMS)  (Shchepetkin  and  McWilliams, 
2005),  and  a  comparison  of  three  simu 
lations  in  April  2004  is  displayed  in 
Figure  7:  WFS  ROMS  nested  in  climato¬ 
logical  temperature  (T)  and  salinity  (S) 
fields,  WFS  ROMS  nested  in  HYCOM 
T  and  S,  and  WFS  ROMS  fully  nested 
in  HYCOM.  April  2004  was  chosen 
because  the  Loop  Current  was  very  stable 
during  this  month  (i.e„  there  was  no 
eddy  shedding  and  the  Loop  Current 
trajectory  was  stable  and  close  to  the 
climatology).  A  priori,  the  model  forced 
with  climatology  should  he  able  to  repro¬ 
duce  the  deep-ocean  currents  reasonably 
well  in  this  situation.  In  the  model  run 
with  climatological  T  and  S  boundary 
conditions,  the  Loop  Current  is  too 
weak.  With  an  SSH  maximum  of  0.3  m 
(upper  left  panel  of  Figure  7),  the  SSH  is 
significantly  lower  than  the  maximum 
derived  from  altimetry  (0.7  m;  lower 
right  panel  of  Figure  7).  By  using 
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HYCOM  T  and  S  (upper  right  panel  of 
Figure  7),  the  model  is  able  to  represent 
the  Loop  Current  more  realistically.  In 
particular,  the  maximum  SSH  (0.6  m) 
is  closer  to  the  altimetry,  but  is  still  too 
low.  By  also  using  HYCOM  velocity 
and  surface  elevation  (lower  left  panel 
of  Figure  7),  the  Loop  Current  strength 
(maximum  SSH  of  .66  m)  comes  closest 
to  the  observations.  This  result  indicates 
that  the  density  field  alone  is  not  suffi¬ 
cient  to  completely  represent  the  Loop 
Current  transport.  The  model  results 
from  the  three  nesting  configurations 
were  compared  with  a  series  of  moorings 
located  on  the  West  Florida  Shelf  (not 
shown).  The  model  forced  with  clima¬ 
tology  is  too  cold  during  summer  and  has 
the  highest  RMS  error  of  all  three  model 
experiments.  During  summer,  the  model 
experiment  forced  with  climatology 
develops  an  unrealistic  southward  coastal 
current.  The  best  temperature  time  series 
is  obtained  with  the  WFS  ROMS  fully 
nested  in  HYCOM,  which  shows  a  more 
realistic  current  variability.  The  WFS 
ROMS  model  nested  in  HYCOM  is  now 
run  on  a  daily  basis  and  the  full  model 
results  can  be  accessed  and  visualized  at 
http://ocgmodl.marine.usf.edu/WFS/. 

The  second  coastal  application 
consists  of  a  coupled  biophysical 
modeling  system  for  the  South  Atlantic 
Bight  and  the  Gulf  of  Mexico.  This 
effort  is  led  by  Ruoying  He  from  North 
Carolina  State  University.  The  regional 
ocean  model  (Figure  8)  consists 
of  ROMS  and  is  nested  within  the 
HYCOM  ocean  prediction  system;  it 
has  5-km  grid  spacing  and  36  terrain¬ 
following  layers  in  the  vertical.  There 
is  no  ocean  data  assimilation,  but  the 
interior  T  and  S  are  relaxed  to  the 
HYCOM  temperature  and  salinity  with 


Figure  8.  Sea  surface  temperature  on  March  1,  2004.  for  the  regional  S  km  ROMS  domain  nested 
within  the  HYCOM  ocean  prediction  system  There  is  no  data  assimilation,  but  the  interior  T  and  S 
are  relaxed  to  the  HYCOM  temperature  and  salinity  with  a  relaxation  time  of  30  days. 


a  30-day  relaxation  time  scale.  This 
system  allowed  for  the  generation  of 
fields  of  biological  variables  from  2004 
to  present  using  the  coupled  model 
hindcasts;  Figure  9  shows  an  example  of 
the  modeled  phytoplankton  distribution 
for  August  27,  2005.  The  ocean  model 
results  were  validated  using  independent 
observations;  Figure  10  gives  an  example 
for  the  University  of  South  Florida’s  WFS 
station  09  near  the  Florida  Keys  when 
hurricane  Katrina  crossed  South  Florida 
in  late  August  2005. 

IMPACT 

In  summary,  25  institutions  (nine  from 
academia,  nine  from  government,  and 
seven  from  the  private  sector)  and 
approximately  60  scientists  worked 
successfully  as  a  team  to  provide  an 
operational  product  that  is  useful  not 


only  for  Navy  and  NOAA  applications 
but  also  for  in  depth  scientific  research 
and  commercial  applications.  This 
effort  would  have  not  been  possible 
without  NOPP  sponsorship.  A  measure 
of  success  is  provided  by  the  interest 
generated  outside  the  partnership.  For 
example,  Kelly  et  al.  (2007)  evaluated  the 
performance  of  a  non-data-assimilative 
Pacific  Ocean  HYCOM  in  a  region 
containing  the  Kuroshio  Extension.  They 
found  that  the  HYCOM  upper  ocean 
heat  budget  is  similar  to  a  diagnostic 
heat  budget  inferred  from  observations 
in  that  the  dominant  contribution  is 
from  lateral  fluxes,  but  that  the  advec- 
tion  fluctuations  are  much  larger  in  the 
model.  Commercially,  the  HYCOM 
nowcasts  and  forecasts  are  often  used  by 
private  companies  to  provide  ocean  state 
information  that  is  critically  important 
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Figuie  9.  Modeled  phytoplankton  distribution  (mmol  m  -)  on  August  27,  200S.  The  pink 
dot  is  the  location  of  Hurricane  Katrina. 
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to  the  cruising,  shipping,  fishing  and 
workboat  fleets.  For  example,  detailed 
surface  current  information  derived  from 
HYCOM  is  summarized  by  OCENS 
(Ocean  and  Coastal  ENviromental 
Sensing,  http://www.ocens.com)  for  their 
customers  who  continually  provide  posi¬ 
tive  feedback  on  their  accuracy. 

OUTLOOK 

The  long-term  goals  of  the  HYCOM 
consortium  for  the  global  domain  are 
to  (a)  add  3-D  and  4-D  variational 
data  assimilation,  (b)  increase  the 
horizontal  resolution  of  the  global 
domain  to  1/25°,  (c)  implement  two-way 
nesting,  (d)  implement  zero  depth 
coastlines  with  wetting  and  drying, 
and  (e)  include  tides.  The  scientific 
goals  include,  but  are  not  be  limited  to: 

(a)  evaluation  of  the  internal  tides  repre¬ 
sentation  in  support  of  field  programs, 

(b)  evaluation  of  the  global  model’s 
ability  to  provide  boundary  condi¬ 
tions  to  very  high  resolution  coastal 
models,  (c)  interaction  of  the  open 
ocean  with  ice,  (d)  shelf-deep  ocean 
interactions,  (e)  upper-ocean  physics 
including  mixed  layer/sonic  depth 
representation,  and  (f)  mixing  processes. 
Other  research  activities  will  focus 

on  coupled  ocean -wave-atmosphere 
prediction;  bio-geo-chemical-optical 
and  tracer/contaminant  prediction; 
ecosystem  analysis  and  prediction;  and 
Earth  system  prediction  (i.e.,  coupled 
atmosphere-ocean-ice-land). 


Figure  10.  Modeled  velocity  fields  versus  observations  at  the  University  of  South  Florida’s  West  Florida 
Shelf  station  C19  during  the  passage  of  Hurricane  Katrina  in  late  August  2005.  The  modeled  tempera¬ 
ture  (lower  panel)  and  observed  surface  winds  (upper  panel)  are  also  displayed. 
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